Dynamical electronic nematicity from Mott physics 
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Very large anisotropies in transport quantities have been observed in the presence of very small 
in-plane structural anisotropy in many strongly correlated electron materials. By studying the two- 
dimensional Hubbard model with dynamical-mean-field theory for clusters, we show that such large 
anisotropies can be induced without static stripe order if the interaction is large enough to yield 
a Mott transition. Anisotropy decreases at large frequency. The maximum effect on conductivity 
anisotropy occurs in the underdoped regime, as observed in high temperature superconductors. 

PACS numbers: 74.20.-z, 71.10.-w, 74.25. Jb, 74.72.-h 



Nematics are translationally invariant but sponta- 
neously break rotational symmetry. They have been 
observed long ago in the field of liquid crystals. In 
the last decade, quantum analogs of these phases have 
been discovered in electron fluids. Quantum-Hall sys- 
tems and many strongly correlated electron materi- 
als with very small structural anisotropy in the plane, 
such as Sr3Ru207, cuprates and pnictides, display very 
large anisotropics^— in properties such as dc (Ref. |j) 
and infrared conductivity^ spin fluctuations, 6 or Nernst 
signal. 7 Microscopic intraunit-cell nematicity in cuprates 
has also been observed^ 

Much attention has been focused on whether nematic 
symmetry breaking in correlated electron materials, such 
as the cuprates, originates from the electron fluid or not. 
In fact, symmetry breaking in the electron fluid necessar- 
ily implies symmetry breaking in the structure and vice 
versa. 2 Recent studies show that a small band-structure 
anisotropy can lead to amplification of the anisotropy in 
low-energy excitations, even when spontaneous nematic 
symmetry breaking does not occur in the electron fluid. 
Slave-boson mean-field theory 9 calculations, for exam- 
ple, suggest that the anisotropy in the magnetic exci- 
tation spectrum 6 - can be enhanced by the proximity to 
a correlation-induced Pomeranchuk instability. Another 
recent paper proposes an enhancement of anisotropy due 
to relaxation of frustration on the Kagome lattice which 
originates from the renormalization of the band structure 
near the correlation-induced metal-insulator transition, 
known as the Mott transition^ Both these studies focus 
on the renormalized shape of the Fermi surface or on the 
electronic dispersion relation. 

Using quantum cluster methods, here we find instead 
that a small orthorhombic distortion of hopping in a 
model for cuprate superconductors leads to a large ne- 
matic anisotropy of the self-energy at small doping and 
large enough interaction strength. The large anisotropy 
that we find also decreases when frequency is too large. 
For all these reasons the name dynamical nematicity is 
appropriate. The self-energy anisotropy dominates over 
the Fermi-surface deformation and appears in the pseu- 
dogap regime 4 - at low temperature as a finite-doping sig- 



nature of the Mott transition^ No additional sponta- 
neously generated one-dimensional structures that could 
also break translational symmetry, such as stripes, are 
needed. This explains why the anisotropy that is easily 
observed in transport quantities is often hard to asso- 
ciate to static stripes (as in YBa2Cu30 y ), even though 
such stripes are observed in some materials X Our results 
are in qualitative agreement with an gle- resolved photoe- 
mission spectroscopy (ARPES) (Ref. Qjl) and conductiv- 
ity measurements 4 in highly underdoped YBa2Cu3 0j,. 
There, large conductivity anisotropies are observed with- 
out major contributions from the CuO chains 4 - and at 
temperatures above the appearance of quasielastic in- 
commensurate peaks^ 

We consider the two-dimensional Hubbard Model as 
the simplest one that exhibits the physics of electronic 
correlations 

h = - ]T Ujdldjv + uJ2 (i) 

ijtr i 

Here, di a is the annihilation operator for an electron with 
spin a at site i and U is the screened repulsive Coulomb 
interaction. The band structure part is described by tij 
which includes the transfer integral for i ^ j and the 
chemical potential [i for i = j. The next- nearest- neighbor 
transfer integral is t' and the band anisotropy is intro- 
duced via the nearest-neighbor transfer integral t along 
the x and y directions as t XtV = t(l ± Sq/2). We ne- 
glect small kinematic effects associated with anisotropy 
in lattice constants and take orthorhombicity into ac- 
count only dynamically through hopping. 

We employ the cellular dy namical-mcan-field theory 
(CDMFT) (Refs. and HBhand the dynamic cluster 
approximation (DCA) fRefsJltil and U3) at zero temper- 
ature T; these methods are capable of capturing the full 
dynamics (i.e., the frequency dependence of the spectral 
function) and the short-ranged spatial correlations be- 
yond the single-site dynamical-mean-field theory. Both 
the CDMFT and DCA map the bulk lattice problem onto 
an effective Anderson model describing a cluster embed- 
ded in a bath of noninteracting electrons. The short- 
ranged dynamical correlations are treated exactly within 
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the cluster. In DCA the embedded cluster is, in general, 
translationally invariant but not in CDMFT. 

In this study, the 2x2 plaquette coupled to bath or- 
bitals is parameterized as 

4 4 
Hdust = - X) E tiAadi" + U J2 d \t d ^4i d H 

+EE + E( V Ki4^Ka + h.c) 1.(2) 



a=l Kcr 



i=l 



The first two terms describe the interacting sites on 
the cluster. For CDMFT, = Uj while for DCA, 
tfj = ^tjj for nearest- neighbor hoppings and t^- — 
■^tij for second-neighbor hoppings. 17 In both cases, 
= /i and the interaction term is the same as in 
Eq. (UJ). The third and the fourth terms represent 
bath-orbital levels and cluster-bath hybridizations, re- 
spectively. Bath orbitals (a) are classified by the clus- 
ter momenta K = (0, 0), (it, 0), (0, tt), (tt, tt) (Refs. [H 
and fl9l ) and the impurity-bath hybridization reflects the 
symmetry as — v^e lK ' ri . The Anderson model, 
Eq. @, is solved using the Lanczos exact diagonaliza- 
tion technique i^£r— This technique limits the number of 
bath orbitals (N a = 2) but allows us to access the dy- 
namical quantities directly on the Matsubara, iuj n , or 
real- frequency axis uj + ir\. 

The Anderson model is subjected to a self-consistency 
condition, relating its one-particle cluster Green's func- 
tion G^ p {iuj n ) = -jdre^ T (T T di(T)d](0)) to the lat- 
tice Green's function of the original model as 



G*"*(«j n ) = N -J{0) G(k,iu n ), (3) 

~ ~ r . ~ -1-1 

G(k,iu n ) = (iu> n + fi) I - t(k) - S(iuj n ) , (4) 



where E(iw n ) = G^T (iui n )- 1 - G %mp (iujn)- 1 . Here, 

= N- 1 Y^ K e*i K+eS )<"- r ^e K . i describes the 
hopping between the clusters covering the original lattice 
with ?? = 1 for CDMFT and d = for DCA,! 7 . k are wave 
vectors in the reduced Brillouin zone, £k = —2(t x cos k x + 
t y cos k y + 2t' cos k x cos k y ) the non-interacting dispersion 
relation of the lattice, and N c = 4 for our 2x2 plaquette. 
The self-consistency Eq. © is solved by iteration, recom- 
puting the bath Green's function G^™f , and thus new 
bath parameters, at each iteration by minimizing with 
a conjugate gradient algorithm a distance function that 
includes frequency dependence on the discrete Matsub- 
ara frequency w„ = (2n — 1)tt/(3, with the small energy 
cut-off (3t = 50, a frequency weighting |l/a; n |, and a large 
frequency cut-off of 10t i 21 ' 23 

The cluster self-energy S (iuj n ) enters a CDMFT lat- 
tice Green's function G(k, iuj n ) with larger unit cell. 
From the self-energy functional point of view2» 4 it is 
natural to focus on diagonal quantities. This corre- 
sponds to Green's function periodization G(k, iuj n ) = 



k/n 
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FIG. 1: (Color online) Spectral function at the Fermi level 
A(k, uj — 0) in the first quadrant of the Brillouin zone ob- 
tained with CDMFT, U = 10, t' = -0.3t. The self-energy is 
computed directly in real frequency with a small imaginary 
part irj with 77 = O.lt. The carrier density is (a) = 0.62, 
(b) 0.81, (c) 0.89 and (d) 0.93. The maximum of the spectral 
intensity A max is given by (a) 0.9, (b) 0.6, (c) 0.4, and (d) 
0.3 in units of 1/t. The thin white lines indicate the location 
of the maxima of ^4(k, ui = 0) that define the Fermi surface. 
The band anisotropy is So = 0.04. 



^c" 1 E 4J e- ,;k '( r '- r ^ G tj (k,iuj n )g. Other periodization 
schemes have been proposed such as cluster cumulant 26 
and self-energy^ 7 - periodizations. In DCA, the conduc- 
tivity can be computed including vertex corrections and 
without restoring the translational lattice invariance of 
the self-energyi2£ We will show that the results we find 
with CDMFT and Green's function periodization are 
qualitatively consistent with DCA. This validates our 
approach. In CDMFT, a general procedure for vertex 
corrections has not been developed yet. 

We use t' = — Q.3t and U = lOt, as parameters that 
are appropriate for cuprates. The U dependence of the 
results for the conductivity anisotropy will be given to 
illustrate the physics. We take as intrinsic orthorhom- 
bic anisotropy in the nearest-neighbor hopping <5o = 0.04 
which is small, like in most cuprates 4^ For So = 0, in 
agreement with previous quantum cluster calculations, 29 
we do not find spontaneous breaking from G4 to Ci sym- 
metry. In other words, we do not find a nematic instabil- 
ity. This kind of spontaneous symmetry breaking is also 
often referred to as the Pomeranchuk instability or effect. 
It is possible that in larger clusters, spontaneous symme- 
try breaking occurs at the van Hove doping, as found at 
weak coupling. 30 Indeed, near this doping, (N = 0.727 
at So = 0) we had difficulty in obtaining well-converged 
solutions when <5o i= 0. 

Figure [1] illustrates the strong anisotropy in j4(k,cij = 
0) = — ilmG(k, w = 0) that emerges as one approaches 
half filling. The results are displayed in the first quad- 
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FIG. 2: (Color online) Angle-resolved spectral function 
A(k, cj) near the Fermi level obtained from CDMFT with the 
same parameters as Fig.l. (a) N = 0.81 and (b) N = 0.93. 
Spectral intensity is given in units of 1/t. The Fermi level 
uj — is indicated by a broken line. In the insets, the vertical 
axis near the Fermi level is blown up. 



rant of the Brillouin zone. The Fermi surface, denned as 
the location of the maxima of A(h, lo = 0), is shown as a 
thin white line in Fig. 1. It is not strongly anisotropic, 
in other words the real part of the self-energy is not 
strongly affected. On the other hand, if we define over- 
doping as N < 0.85 and undcrdoping as N > 0.85, 22 
then in the latter regime the spectral intensity shows 
strong directional anisotropy; at N = 0.93, for exam- 
ple, the spectral weight near (it, 0) is suppressed by as 
much as 54% compared with that near (0,7r). In other 
words, the pseudogap that normally appears symmetri- 
cally near (7r, 0) and (0, n) is now much more strongly 
anisotropic than the 4% band anisotropy. This reflects a 
strong anisotropy in the self-energy. Recent ARPES ex- 
periments on YBa2Cu30j, have shown similar anisotropy 
in the lightly doped regime*^ 

Let us now turn to the finite frequency regime. In 
Fig. [3J we present the angle-resolved spectral function 
A(k,w) for the overdoped (a) and the underdoped (b) 
regimes. The blow up of the region near the Fermi level in 
the insets shows that the anisotropy is most pronounced 
in the underdoped regime, N = 0.93. The anisotropy 
also extends to progressively higher frequencies as the 
system is underdoped. The strong anisotropy essentially 
disappears at high frequency, reminiscent of recent obser- 
vations on the spin-fluctuation spectrum^ This indicates 
that the nematicity here is a dynamical phenomenon as- 
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FIG. 3: (Color online) Anisotropy in the CDMFT conduc- 
tivity 8a = 2 [(7a; (0) — cr y (0)] / [<Ta;(0) + (T y (0)] as a function of 
filling TV for various values of U and n = O.lt, do = 0.04. 



sociated with the electronic response near the Mott tran- 
sition as N — )• 1. 

The link with Mott physics is seen most clearly in the 
anisotropy of the dc conductivity cr x (y) that we consider 
now. For the CDMFT results, we use the Kubo formula 
with the pcriodized Green's function to compute 
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[lmG(k, 0) 



This formula neglects vertex corrections. 



(5) 

The con- 



ductivities u x (y) extrapolate to infinity like I/77 as the 
small imaginary part 77 tends to zero, as expected in 
a pure metal (Drude peak) at T = 0. We have 
checked however that a y /a x reaches a limit. Never- 
theless, we use 77 = O.lt in Fig. to avoid uncertain- 
ties in extrapolating the conductivity anisotropy S a = 
2 [cr x (0) - a y (0)] I [a x {0) + ov(0)]. With smaller 77 the 
anisotropy is generally larger. 

As a function of filling N and for all U, there is a 
peak around N ~ 0.8, not far from the non-interacting 
5q = van Hove singularity at = 0.727. This is a 
band-structure effect since the Fermi surface opens up as 
if the system had a tendency to be quasi-one-dimensional 
[see Fig. Q] (b)]. We found little 77 dependence of 5 a near 
N = 0.8 because of the large imaginary part of the cluster 
self-energy for that filling^ 

The most interesting results occur close to half fill- 
ing. For values of U below the critical U& for the Mott 
transition^ the anisotropy takes small values consistent 
with the small orthorhombic distortion. At U = 6rj, the 
anisotropy in conductivity is largest and not monotonic 
close to N = 1. This U is just slightly above the critical 
Ud ~ 5.25i for the Mott transition obtained for CDMFT 
with exact diagonalization and t' — 0.— The abrupt 
and large increase of the anisotropy at a finite doping 
close to half filling for all larger values of U is similar 
to the experiment. 4 The CDMFT results are consistent 
with the existence of a first-order transition between two 
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FIG. 4: (Color online) Anisotropy in the DCA conductivity 
for So = 0.04 as a function of filling N for various values of U 
and 7] = O.lt. Open symbols do not include vertex corrections, 
filled symbols do. 

kinds of metal found recently for U larger than U C 2 at 
finite dopingjii The metal closest to half filling displays 
a pseudogap and we have shown in Figs. [1] and [2] that in 
this metallic phase it is extremely sensitive to small or- 
thorhombic distortions. The anisotropy in occupancy of 
K orbitals, at most 12%, does not suffice to explain the 
anisotropy in conductivity. The imaginary part of the 
self-energy plays a prominent role since replacing it by a 
constant in the optical conductivity formula removes the 
large anisotropy. The fact that the enhanced anisotropy 
occurs far from the van Hove singularity suggests that it 
is not simply a weak-coupling effect like the Pomeranchuk 
instability^ 

The above results are confirmed by Fig. [4j ob- 



tained from DCA calculations that include vertex 
corrections i 28 ' 33 ' 34 The value of U c is about 4.5i for 
DCA at t' = 0i2S It is above this value that the large 
anisotropics are present. The critical filling where the 
anisotropy increases does not depend strongly on <5rj. 
Anisotropy scales like Sq at small Sq and is linked with 
the pseudogap so it should appear at comparable Ti^I 
Other interesting physics has been found in DCA at 
U > [/ e ) 35 i 38 including a critical filling similar to ours 
at U = 8t beyond which a local bond-order susceptibility 
diverges at T — >■ 022- 

We have shown that, in the presence of interactions 
larger than the critical U for the Mott transition, small 
orthorhombicity can lead to very large anisotropies in the 
pseudogap and conductivity in the underdoped regime 
close to half filling without further symmetry break- 
ing. Transport anisotropy should be much more sensi- 
tive to uniaxial pressure in the underdoped than in the 
overdoped regime. Mott-induced dynamical nematicity 
should influence all other transport properties. The com- 
petition between this effect and other symmetry-broken 
phases would be worth investigating. Preliminary study 
reveals that, in the presence of small band anisotropy, 
the superconducting order parameter becomes highly 
anisotropic, A x ^ — reminiscent of experimenti 41 - 
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